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Abstract 

The local Casimir energy density for a massless scalar field associated with step-function poten- 
tials in a 3 + 1 dimensional spherical geometry is considered. The potential is chosen to be zero 
except in a shell of thickness 6, where it has height h, with the constraint h5 = 1. In the limit 
of zero thickness, an ideal <5-function shell is recovered. In this limit, the behavior of the energy 
density as the surface of the shell is approached is studied in both the strong and weak coupling 
regimes. The former case corresponds to the well-known Dirichlet shell limit. New results, which 
shed light on the nature of surface divergences and on the energy contained within the shell, are 
obtained in the weak coupling limit, and for a shell of finite thickness. In the case of zero thick- 
ness, the energy has a contribution not only from the local energy density, but from an energy 
term residing entirely on the surface. It is shown that the latter coincides with the integrated local 
energy density within the shell. We also study the dependence of local and global quantities on 
the conformal parameter. In particular new insight is provided on the reason for the divergence in 
00 ■ the global Casimir energy in third order in the coupling. 
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I. INTRODUCTION 



The subject of local energy density associated with the confinement of quantum fields 
by surfaces has a rather long history. For example, Brown and Maclay p| computed the 
vacuum expectation value of the electromagnetic energy-momentum tensor between two 
parallel perfectly conducting plates, which is twice that of a conformally coupled scalar field 
satisfying Dirichlet or Neumann boundary conditions on the plates separated by a distance 
a in the x direction, namely 

( T H = Ti ^diag (-1,-3, 1,1), (1.1) 

which corresponds precisely to the attractive energy or pressure found by Casimir in the 
same situation |2j . If a nonconformal scalar stress tensor is used, a position-dependent term 
in the stress tensor appears, which does not contribute to either the total energy or the 
pressure on the plates 0, 0| • 

Local surface divergences were first discussed for arbitrary smooth boundaries by Deutsch 
and Candelas 0|. They found cubic divergences in the energy density as one approaches 
the surface; for example, outside a Dirichlet sphere (that is, for a conformally-coupled scalar 
field satisfying Dirichlet boundary conditions on the surface) the energy density diverges as 

11 . . 

r — > a+ : u ~ - — — , (1-2) 

3607r 2 a(r-a) 3 ' v ; 

where a is the radius of the sphere, and r is the distance from the center of the sphere. This 
raises the question: How can it be that the total Casimir energy of a Dirichlet sphere (or a 
perfectly conducting sphere in electrodynamics) is finite? The electromagnetic case is the 
well-known one first calculated by Boyer [f| 

£E „ = (L3) 
a 

while the scalar case was first worked out by Bender and Milton 0], 

fS= O£02817 
a 

Thus there has been a feeling since the time of Deutsch and Candelas that there was some- 
thing suspect about the calculations of Casimir self energies of ideal closed boundaries. (We 
note that there is now a proof that any such smooth perfectly conducting boundary pos- 
sesses a finite electromagnetic Casimir energy 0|. Whether such an idealized limit is physical 
is, of course, another question.) This suspicion has been recently intensified by a series of 
talks and papers by Graham et al. |9| . The essential outcome of their analysis is that for a 
5-function sphere, described by the following Lagrangian for a massless scalar field, 



C = —d^dy - ~5(r - a)0 2 , (1.5) 



a divergence occurs in third order in A. (They claimed a divergence in second order, but that 
was spurious BE EH-) This diver gence in fact was discovered much earlier by Bordag, 
Kirsten, and Vassilevich [l2|], and possible ways of dealing with it have been suggested 
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[HI Hl|. (Recently, the effects of the spherical 5-function potential at finite temperature 
have been considered by McKenzie-Smith and Naylor 0|.) Objections complementary to 
those of Ref. have also been voiced by Barton ^6[, all of which raise doubts concerning 
the physical relevance of results such as those in Eqs. (|1.3|) and (jl.4|) . 

Clearly, then, there are issues here unresolved. The purpose of the present paper is to 
examine the local energy density for a scalar field in the presence of a spherically symmetric 
background, but as suggested by Ref. [^| not so singular as a delta function. So we consider a 
spherical shell, with a finite step potential. In the limit as the height goes to infinity and the 
width to zero we recover the 5-function sphere. This enables us to explore how the quantum 
vacuum energy of the shell itself contributes to the total energy of the configuration. In the 
next section, we will construct the Green's function for such a sphere, and in the following 
section the local energy density (for an arbitrarily coupled scalar) inside, outside, and within 
the region of the shell. In Sec. I1VI we will consider the thin shell limit and examine how 
the energy density diverges as the surface is approached. As expected, the divergence in 
the energy density as the surface is approached is weakest for the conformally coupled case. 
We also study how the planar divergences within and outside a thick shell pass over to the 
divergences exterior to a thin shell. In Sec. IS we will compute the energy of each region, 
and thereby determine the total energy of the spherical geometry. There we will see that the 
former is constructed not only from the local energy density, which diverges as the surface is 
approached, but also contains an energy lying entirely on the surfaces bounding the region, 
except for a special, nonconformal value of the conformal parameter. For finite thickness 
of the shell, the surface energies cancel across each surface in constructing the total energy. 
However, the thin shell limit corresponds to a singular potential, for which there is a net 
effective surface energy for the sum of the interior and exterior regions of the shell which is 
identified as the integrated energy density within the shell. 



II. GREEN'S FUNCTION FOR A SPHERE 



As discussed in Refs. |10|, [ll|, [17( for parallel planes, we consider the potential 

Ant = -^0V(r), (2.1a) 

where 

0, r < a_ 

a(r) = ■{ h, a_ < r < a+, (2.1b) 
0, a + < r. 



- > 



Here a± = a ± 5/2, and we set hS = 1. Compared to Refs. jlOL we have changed the 
dimensions of A so that the total energy of interaction does not explicitly refer to the radius 
a. In the limit as 5 —>■ (or h — > oo) we recover the 5-function sphere considered first in 
Ref. Hi. 



A straightforward solution of the Green's function equation, for a massless particle, with 



K 2 = —LU 2 , 



(-V 2 + k 2 + ±A Q(r, r') = 5(r - r'), (2.2) 
is given in terms of the reduced Green's function, 

£(r,r') = ^^y)*UM)^n(^0O- (2-3) 



tin 



3 



The latter is expressed in terms of the modified Ricatti-Bessel functions, 



= J^I l+ i/ 2 (x), ei(x) = )J^-K l+1/2 (x) 



as follows, outside of the shell, 



r, r < a_ : g x 



r.r > Oj 



9l 



1 

Krr' 



Sj(«r<)ej(K;r>) - ^ei(K,r)ei(nr') 



(2.4) 

(2.5a) 
(2.5b) 



where the denominator is 

S = [KS;(/ta_)e;(K'a_) — K,'si(Ka-)e' l (K , a-)][K'ei(Ka+)s' l (K'a+) — «;e;(fi;a + )s;(fi;'a + )] 

— [Ks'^KaJjsi^' a) — k' si(na)s\(K! aJ)][K! ei(na + )e[(K' a + ) — /cej(/ca+)ej(«'a-i-)]. 

(2.6) 

The numerator in (|2.5aj) . H, is obtained from H by replacing sj(/ca_) — ► ej(/ca_): 

5 = [KeJ(Ka_)e;(K'a_) — /?'e;(fi;a_)e[(fi;'a_)][K / ei(fi:a + )s;(K / a + ) — /«eJ(fi;a + )si(K / a + )] 

— [«;eJ(fi;a_)s/(K / a_) — K , e/(fi;a_)sJ(fi: / a_)][K / e/(fi:a + )eJ(fi: / a + ) — «;ej(«a + )ej(« / a + )], 

(2.7) 

while that in (|2.5b|) . S, is obtained from H by replacing ei(na+) — > sj(/ca+): 

H = [/€sj(fi;a_)ei(fi; / a_) — K , si(Ka-)e' l (K , a-)][K, f si(Ka+)s f i(K, f a+) — KS;(/ta + )s/(/t'a + )] 

— [/csj(/ca_)sj(/c'a_) — /t / s/(/ta_)sJ(/t / a_)][ft / s;(fta + )e;(K / a + ) — /tS;(/ta + )e/(/t'a + )]. 

(2.8) 



Here k' = y k 2 + \h/ a 2 . It might be noticed that under the interchange 

si(na T ) <-> ej(/ca ± ), 



(2.9) 



and similarly for functions depending on k', 5 is unchanged, while H <-> H. 

The Green's function within the shell has a somewhat more complicated form. For 



a — | < r, r' < a + ~ we find 



K'rr 



7 — Sz(K / r < )e/(K / r > ) - - 



A[s/(K / r)e«(fi;V / ) + sj(/c / r / )e i («V)] 



+ Bsi(K'r')si(n'r) + Cei{K'r')ei{K'r) 



(2.10) 



where 



A = [fi;e;(fi;a + )ej(K / a + ) — K / e/(/€a + )e;(fi; / a + )][KS;(fi;a_)s;(K / a_) — /t's;(Ka_)sJ(/t'a„)], 

(2.11a) 

B = — [/s;e|(fi;a + )e/(/« / a + ) — /t'e/(fi;a + )e;(fi; / a + )] [KsJ(Ka_)e;(fi;'a_) — K/sj(/ca_)eJ(«'a_)], 

(2.11b) 

C = — [/te;(fi;a + )si(fi; / a + ) — K / ez(Ka + )sJ(ft / a + )] [KsJ(/ta_)sz(/t'a_) — Ac'sj(Ka_)sJ(Ac'a_)]. 

(2.11c) 
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A check of this occurs when we take the interior radius to zero, a_ 



0, for then we recover 



the known result for (the TE part of) a solid dielectric ball Then si(na_) — > 0, and 

the A and C coefficients above vanish. The denominator also simplifies, and we obtain the 
following form for the Green's function within and without the ball, of radius a + now called 
a: 



where 



gi(r, r') 



K'rr 
1 



7[sj(K / r<)e l (/c'r>) - &si(K , r)s l (K'r 1 )], r,r' < a, 



nrr 



7 [s / (/tr < )e / (fi;r > ) - fiej(«r)ej(Kr')], r,r' > 



n 



n 



lim 



B 



lim 





[na 


)e 


'(n'a) 


— Ke 


{na 


)ej 


[n'a) 


K'ei 


' Ka 


)s 


(n'a) 


— ne 


(na 


)si 


[n'a) 


k'si( 


na] 


s\ 


(k'o) 


— Ks\ 


(na] 


Sl( 


n'a) 


n'ei( 


na] 


s\ 


(n'a) 


— ne\ 


(na] 


si( 


k'cl) 



(2.12a) 
(2.12b) 

(2.13a) 
(2.13b) 



The required symmetry between the inside and outside forms of the energy, obtained by 
interchanging si «-> e\ and k <-» k', is now manifest. This result for the Green's function 
is exactly that found in Ref. ^| for the TE (scalar) part of the dielectric Green's function 
with fi= 1, k! = ■sJe! Jen. 

We will now use the Green's function to compute the energy density. 



III. ENERGY DENSITY 

We can calculate the local energy density from the stress tensor 

= d"^d v (j) - g^C - f - g^d 2 )<f) 2 
from which the energy density follows: 

A 



T 



oo 



& 



OioO. 



Vd>- Vc 



2,2 



(3.1) 



(3.2) 



where the conformal value is given by £ = 1/6. To obtain the vacuum expectation values at 
one loop, we use the identification with the classical Green's function 



{<j>{x)<j>{x')) = -G(x,x'). 



The energy density thus is, within or outside the shell, 



<T 00 ) 



2i 



,0o/0 



d"d' 



V V 



G(x,x') 



V 2 G(x,x). 



(3.3) 



(3.4) 



When we insert the partial wave decomposition of the Green's function ([2.3)1 . we en- 
counter, in terms of the angle 7 between the two directions (6, ip), (6', ip'), 



cos 7 = cos 9 cos 6' + sin 6 sin 6' cos(y? — <p') . 



(3.5) 
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the evaluation 



V • VP; (cos 7) 



9',<j>=4>' 



>?<D = ^ 



(3.6) 



Therefore, the expression for the energy density is immediately reduced to (inside or outside 
the shell, but not within it) 



(T 00 > 



dn A 21 + 1 



E 



2tt ^ 4tt 

z=o 



-k 2 + d r d r i + 



Z(Z + 1) 



2£——r 2 —g (r r)\ 
r 2 dr dr 'J 



Note that there has been no need of some dubious argument (such as appears in Ref. jlll |. 
Eq. (4.21)) concerning partial integration in the angular coordinates. 

We insert the Green's function in the exterior region, but delete the free part, 

1 



9i 



nrr 



7 s i (/tr < )e i (/tr > ), 



(3.8) 



which corresponds to the bulk energy which would be present if either medium filled all of 
space, leaving us with for r > a + 

21 + IS f e 2 (nr) 



u(r) 



-(1-40 



dn 



E 



2lX ^ 47T 

1=0 



-K 



; i+4e | | 1 

ly'l tyl 



(3.9) 



- — eiinry^Kr) + —e[ («r) 

Inside the shell, r < a_, the energy is given by a similar expression obtained from Eq. ()3.9|) 
by replacing ei by si and S by S. 

Within the shell, a_ < r < a + , the energy density is given by a somewhat more compli- 
cated formula: 

/2 l + 4£ l(l + l) + l 2\h 1 



n(r) 



40 



G?K ^ 2Z + 1 1 1 

27T ^ 47T K'r 2 S 
Z=0 

2/, 



l-4£ 



+ 



a 2 l-4£ 



x[2As/(K/r)ej(KV) + Bsf^'r) + Cef(K'r)} 
2k' 



— [A(si(n' V)ej(/t' V) + e/(K/r)sj(/«V)) + 5s/(/tV)sJ(/tV) + Cez(fi/r)eJ(fi;V)] 



+ K ,2 [2A S ;(fi: , r)e;(fi:V) + Ssf («V) + Cef («'r)] 
Here A, B, and C are the coefficients given in Eq. ([2. lip . 



(3.10) 



IV. SURFACE DIVERGENCES IN THE ENERGY DENSITY 
A. Exterior divergences for S = 

We want to examine the singularity structure as r — > a+ from the outside. For this 
purpose we use the leading uniform asymptotic expansion, I — > 00, 

e,(x) ~ v^te-^, s,(x)~i>/ite^, 

ej(x) )=e~^ 4( x )~\-±=<f, (4.1) 
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where [y = I + 1/2) 



x = uz, t = (l + . 2 )- 1/2 , d J = - + - (4-2) 

dz zt 

Let us first consider the thin shell limit, 8 — > 0, hS — 1, where it is easy to check that 

Ka n ! (4.3) 



which is exactly the coefficient occurring in the ^-function potential 1)1.5)1 . There are two 
simple limits of this, strong and weak coupling: 

A 5 sAna) 

- oo : - -> -L (, (4.4a) 
a c e/(Avaj 

- - : 1 - -A^M, (4.4b) 
a n Aca 

if we assume that the relevant scale of k is 1/a, since we expect that the significant values 
of k are determined by the argument of the Bessel functions in Eq. ()3.9)1 . 

In either case, we carry out the asymptotic sum over angular momentum using Eq. ()4.1)l 
and 

oo 1 

y e -^ = — — . (4.5) 

2sinh| v ' 

z=n 2 



Here (r 



X 



r 



drj r — a 2r — a 



T)(z)—T)[z-) m2z—(z) = . (4.6) 



dz r t r 

The remaining integrals over z are elementary, and in this way we find that the leading 
divergences in Eq. (J3.9)) are as r — > a+, 

X 1 l-6f lAn , 

> oo : it ~ — , (4.7a) 

a (r — a) 4 

--0: uW ~ f-^V r ^ 4 -") (!- 6Q f-g-V"", n<4, (4.7b) 
a \ ay 967Ha 4 \r — a J 

where the latter is the leading divergence in order n. These results clearly seem to demon- 
strate the virtue of the conformal value of £ = 1/6; but see below. (The value for the Dirichlet 
sphere f)4.7a)l first appeared in Ref . |Bj ; it recently was rederived in Ref . [l^| , where, however, 
the subdominant term, the leading term if £ = 1/6, namely 1)4.12)1 . was not calculated. Of 
course, this result is the same as the surface divergence encountered for parallel Dirichlet 
plates 

Thus, for £ = 1/6 we must keep subleading terms. This includes keeping the subdominant 
term in x, 

2 r — a I r — a 



-M — J > ( 4 - 8 ) 
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the distinction between t(z) and t = t(z = za/r), 

n r — a , 

ztxizt-t 3 z , (4.9) 

r 

as well as the next term in the uniform asymptotic expansion of the Bessel functions, 

Sl (x) ~ -Vzie"" (l + + 0(u- 2 )^j , (4.10a) 

ei{x) ~ Vzie'^ (l - ^u x {t) + 0(z/- 2 )^) , (4.10b) 

s\{x) ~ i-L e ^ + + (^ 2 )) , (4.10c) 

e{(s) ~ e -»i ^ _ + 0(z7 -^ > (4 . 10 d) 

where 

, , 3t - 5t 3 , , 3t + 7t 3 , , , , , 

Including all this, it is straightforward to recover the well-known result (jl.2|) [Hf for strong 
coupling (Dirichlet boundary conditions): 

- -»• oo : m~ — 1 - . 1 — — , (4.12) 
a 3607H ayr — a) 6 

Following the same process for weak coupling, we find that the leading divergence in order 
n, 1 < n < 3, is (r — ► a±) 

/ \\ n i i 

A^O: ^ — — -— ^— ; r=— (n-l)(n + 2)r(3-n). (4.13) 

Va 2 / 14407T 2 a(a - r) 3 ~ n v /v 7 v 7 v 7 

Note that the subleading 0(A) term again vanishes. Both Eqs. (|4.12j) and (J4.13j) apply for 
the conformal value £ = 1/6. 



B. Divergences within and outside the shell for 5^0 

Now we take finite 5, and examine the energy density (J3.10)) within the shell. Again, 
the singularities in the energy density as r approaches a + from within the shell are revealed 
by examining the behavior as k, k' — > oo. From the leading uniform asymptotic expansion 
(|4.ip . the denominator (|2.6|) becomes (I — > oo) 




(4.14) 



S 



where Z- = Kd-/v, z' + = K'a + /u, r]' + = T](z' + ), etc. Because r)' + > r)'_, it is clear that the first 
term in Eq. (|4.14|) dominates, so as I — > oo, 



E ~ e^-^'^-K 2 z ^~ z ^ ( i + !z ) ( i + !+ ) . (4.15) 



4 y z_t-z + t+ \ t'_j v 

As for the numerator in Eq. 1)3.10)1 it is clear from Eq. ()4.1)) that the terms proportional to 
sf(n'r) and sf(n'r) dominate, and so again we need only examine the B coefficient 1)2.11 b)) . 
which is approximated by 



fl ~e^W^^(l + y(l-£). (4,6) 



Thus the leading coefficient in the Green's function or the energy is 



B A-f 

2e -W + _J± (4.17) 



1 4- - 

1 + t' 



+ 



Now we must relate t' + and £ + . t + — (1 + ^i) 1 ' 2 , and 

<< = 1 



,/2„2 



1 / Xh a 2 . 



l + zi 



XhaA 



2v 2 a 2 + ' 



tl , (4.18) 



again as v — > oo. Thus, in this approximation the energy density f)3. lUj) in the shell is 

Xh ai 1 f°° , 



where we have used the evaluation (J4.5j) . or 

oo v v 2 

S>~~fer-)- <«<» 



. 2 r — a+ , 

Then in the limit as r approaches the outer surface 

-< r °°>- 9 6 ^K->f (1 - 60 - r - a +-- (4 - 21) 

Exactly the same result (with a + replaced by a J) is obtained as r approaches the inner 
surface of the shell from within the shell, r — > a_+. 



That Eq. (j4.21j) is exactly the expected result is seen by recalling from Ref. [lOj, [ll|, [l 



the divergence in the local energy density encountered as one approaches the surface of a 
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slab, centered on the origin, of height h and thickness 8 with two transverse dimensions from 
the inside, Eq. (2.52) of Ref. [uj (or Eq. (91) of Ref. 0): (Xh/a 2 -> Xh there) 

This was obtained, as here, by looking at the most divergent part of the k integral, and 
using k 2 ^ Xh. On physical grounds, we must encounter the same divergence here, because 
if we are very close to the surface locally it looks flat. 

We obtain the opposite behavior outside the outer surface, where the same asymptotic 
analysis implies 

E „ -.l^—f. „ i e 2^ + ^^4 (4.23) 
S 2 l + it 2 2v 2 2a 2 K J 

l + 

This implies exactly the negative of the divergence found in Eq. (j4.21j) as r approaches the 
surface from the outside. (Again, the change of sign was found outside the slab [lfl H^|.) 

This, obviously, is quite different from the surface divergence we encounter in Eq. (|4.7bJ) — 
it has a different dependence on the distance from the surface, and furthermore it depends 
explicitly on h. Unlike the curvature divergences found in the previous subsection, the 
divergences found here are independent of curvature [Xh/a 2 is simply the coupling strength 
of the potential in Eq. (j2.1|) ]. it cancels between inside and outside the boundary, and it 
vanishes if the conformal coupling is used. As in the parallel plate situation, it is without 
physical consequence therefore, and we will disregard it in the following, where we wish to 
concentrate on the thin shell, <5 — ^ 0, limit. 

It is, however, easy to provide a heuristic connection between the divergence structure 
seen in Eq. (|4.21j) and that obtained outside a thin shell, Eq. (|4.7bJ) for n = 1. The latter 
is obtained from the former by considering the sum of the energy densities due to the inner 
and outer surfaces, a± ^> r — a± ^> 8, 



u(r -»• a++) + u(r -»• a_+) = — ^ 2 (1 - 6£) 



1 1 

+ 



(r — a_) 2 (r — a+) 



A(| - ( *> a + -a_=5^0, (4.24) 



A8n 2 a 2 (r - a) 3 ' 



again using h8 — 1. This is exactly the expected divergence, which therefore has nothing to 
do with curvature, since the same result may be obtained either by this argument, or direct 
calculation, for a thin (5-function) plane slab. 



V. TOTAL ENERGY FOR A THIN SPHERICAL SHELL 
A. Cancellation of divergences 

The results ()4.13|) from Sec. II V Al for the conformally coupled scalar for a thin shell show 
that the inverse linear divergences in the local energy density which occur in either order 
A or A 2 cancel between inside and outside, when one computes the total energy, while the 
divergence encountered at n = 3: 
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is correctly given by 



A 3 r — a 
m 



144vr 2 a 7 a 



as shown by explicit calculation, because 

oo 1 



(5.2) 



(5.3) 



The integral of u^ 3 \ however, is finite, so this does not signal any difficulty with the total 
energy in order A 3 . 



B. Surface energy 



However, as discussed first by Dowker, Kennedy and Critchley |2 



2 11 1 , a nd later elabo- 



rated by Saharian and Romeo 22], and put in a broader context by Fulling [23j], for situations 
when other than Neumann or Dirichlet boundary conditions apply, an additional term must 
be supplied in calculating the energy, a term which resides entirely on the surface. For the 
case of the general stress tensor (j3.1|) . that extra term is 



2i 



dS ■ VG(x, x) 



(5.4) 



where the direction of the normal is out of the region in question, which arises from the T' 
component of the stress tensor, and from d^T^ = 0. The total energy in a given region V 
bounded by a surface 5* is not, therefore, just the integral of the local energy density, but 
has this extra contribution [111 ]: 



£ = /(*)<n + e = iX(*)/^ ( 



r, r e 



(5.5) 



which is independent of £. The latter expression has a rather evident interpretation in terms 
of summing zero-point energies. (For example, see Appendix A of Ref. |24j.) As there, we 
have inserted a temporal point-splitting regulator r as first introduced in Ref. [25| . The limit 
t — > must be taken at the end of the calculation. Of course, the total energy in all space 
is just the sum of the integrated local energy densities in each region, because the sum of 
the inside and outside surface energy for each surface is zero, owing to the continuity of the 
Green's function and its normal derivative across each surface for the nonsingular potential 
(|2.1|) . This is not the case for the singular potential (jl.5|) . 



C. Shell energy 

In the limit of h — > oo for the region in the shell, a_ < r,r' < a + , Eq. (|2.1(Jj) becomes 
[Xh ^> (/ta) 2 , which is the opposite limit to that taken in Sec. II V Bj 



9i 



eiiKdjsAKa) 



2nrr' I + -^ei(Ka)si{Ka) 



cosh 



r — r 



cosh 



r + r' — a + — a_ 



(5.6) 
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(Recall we are disregarding the physically spurious planar surface divergences, discussed in 
Sec. IIVB1 which in any case disappear for the conformal value £ = 1/6.) In the thin shell 
limit this leads to an energy density in the shell nearly independent of r (5 — > 0), leading to 
the energy (y = \x\, e = r/a, r now being the Euclidean time) 

e ° - h<r°> -^-^ f> + 4 £ <twir (57) 

However, we have to include the surface term in the shell, in the thin shell limit, which from 
Eq. JOD is 

r'=r, r=a_|_ 

(5.8) 

r'=r, r=a_ 

where the change in sign in the normal is now encompassed in the limits. (E s exactly cancels 
E s in the thin shell limit: £ s = E s + <£ s — > 0, because the total energy of the shell is given 
by Eq. (fo~5j) 

£ s = \.\ (dr) /^2^(r,r)e— , (5.9) 

2« ./shell J 27r 

which clearly vanishes as the thickness of the shell <5 — 0, because the integral is proportional 
to the volume of the shell. However, we shall see shortly that E s in Eq. ()5.7|) plays a special 
role. [The reader should note that if we set r = expressions such as Eq. (|5.7J1 fail to 
converge. As we shall see shortly, if this point coincident limit is expanded in powers of A, 
divergences in the total energy occur in orders A and A 3 .] 



l-4£ 

2i 



dS 



VG(r,t;r',t 



r) 



D. Total energy of A sphere 

Likewise, if one integrates the interior and exterior energy density in the thin shell limit, 
one gets, for arbitrary £, 

E in + E ont = £ - <E, (5.10) 



where the total energy is that found in Ref. |Uj,ll2|, (y = \x\) 



1 °° 1 f°° d 

£ = > (21 + 1)- / dxy— In 



-I v {y)K v {y) 



(5.11) 



and where the extra term <E is precisely the sum of the surface terms for the inside and 
outside regions. Again, Eq. ()5.10|) is a restatement of Eq. (|5.5jl . and £ is exactly that 
obtained from the integral of the Green's function, as in Eq. ()5.5j) . 

The above can be verified by computing the pressure. This is obtained by evaluating the 
radial-radial component of the stress tensor (J3.1)) , or 



T 

1 v 



2 ]-m 



o 2 - Vi)0 2 , 



(5.12) 
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or more precisely, the discontinuity of the vacuum expectation value of this across the surface: 



P 



(T rr j ; 
1 



(T rr ) c 

oo 

2tt ^ 



1=0 



An 



d r O r i 



K 



1(1+1) 



g(r, r') 



(5.13) 



-r=a+ 



with k = \(\. We note that the £ dependence drops out immediately, and the terms not 
involving radial derivatives of Bessel functions cancel between inside and outside, leaving us 
with, in the thin-shell limit, 



A 1 



vr 2 a 5 2 



oo 



dx) (21 + 1 



1=0 



] y [I v (y)K u (y)]' - h(y)KM 
1 + ±JMK»(y) 



(5.14) 



which, when multiplied by the area of the sphere to give the stress on the surface, is exactly 
that obtained by differentiating the total energy ()5.11|) with respect to the radius (with A 
fixed) , 

df 

~ = 47ra 2 P. (5.15) 
da 

(The scaling law for A was incorrectly stated in Refs. [ull lll| . as was the sign of the pressure.) 

However, there is more to say here. As noted above, the integral of the local energy, inside 
and outside the sphere, is finite perturbatively, because of cancellations between inside and 
outside, for the conformally coupled scalar. But it is well known that divergences occur in 



the total energy at order A 3 BEiQ. These evidently must arise from the surface term 
(15. 4|) . So let us consider the latter, which is given in the outside region by 

d( d 



a 2 (l 



1=0 J 1 



2n dr 



9l(r, r') 



(5.16) 



In the strong coupling limit, there is, of course, no surface term. This is because then 



r,r > a 



9i(r, r') 



1 



Krr 



Si{Kr < )ei{nr >J 



edm) 



ei(nr)ei(Kr') 



(5.17) 



which vanishes on the surface, and has a derivative proportional to the Wronskian. 

In general, in the thin-shell limit, the sum of the inside and outside surface terms is given 

by 



A 



2vra 2 



1 POO _ _ 

;i-4£)-y d*]T(2Z + i) r 



h(y)K v (y) 



(5.18) 



because the free term in the Green's function cancels between inside and out. It is apparent 
that the strong coupling limit of this vanishes, being just a delta function. Perhaps not 
remarkably, this is precisely the same as the integrated local shell energy E s , ()5.7|) . Thus 
the surface energies within and outside the shell regions cancel. (As noted above, this is a 
general result, because the Green's function and its normal derivative are continuous across 
each surface at a + , a_. We have explicitly verified this using (|2.5|) . (|2.1U|) . and the Wronskian 
for the Bessel functions.) For weak coupling, we expand this in powers of A. Again, the 
easiest way to isolate the asymptotic behavior is to use the uniform asymptotic expansion, 



oo : I v (x)KJx) 



t 

2~v' 



(5.19) 
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This yields the following expression for the nth term in the total surface energy, if we set 
e = and regard continuous variable, analytically continuing from Ren > 3, 



e "" ~ -27TJ 1 - 4 «> U) T(tr (2 "" 2 " 1)<(B " 2) - (5 ' 20) 

Note that this expression vanishes for n = 2; we must keep subleading corrections to see the 
order A 2 term in the energy arising from the surface energy. However, for n = 3 we obtain 
for the conformal value, £ = 1/6, 

S<3 ' ~ 5S>>- < 5 - 21 > 



precisely the divergent term in the energy given in Ref. first found by the heat kernel 
calculation of Bordag, Kirsten, and Vassilevich Alternatively if we keep the point-split 
regulator, we find the n = 2 term still vanishes, 

M A 2 i [°° dz 1 , 

<£ (2) ~ / = 0, (5.22) 

24ea 3 ej^ z z 2 + 1 ' 1 ' 

which we argue is imaginary and the integrand odd, if we interpret the integral as its principal 
value. The same phenomenon resulted in the finiteness of the energy for the strong coupling 
limit — the Dirichlet sphere — see Ref. HEl3. The divergence at n = 3 persists: 

£(3) ^ ^Lhie. (5.23) 
127ra 4 v ; 



Thus, as shown in Ref. 11 1, the total energy term 8^ from Eq. (|5.11|) . by expanding the 
logarithm to order A 2 and applying the derivative operator, is unambiguously finite: 

£(2) = (5.24) 
On the other hand, £^ is unambiguously divergent. 



E. Other values of £ 

We have shown, as expected, that for the global energy any value of £ can be used. Fulling 
(23l | has suggested that a value different from the conformal value may have advantages. 
Thus, because the surface (or shell) energy is proportional to 1 — 4£, it vanishes for £ = 1/4, 
as do most of the terms in the energy density (J3.9|) . In that case, however, the surface 
divergences in the local energy density are intensified, which does not seem to account 
for the divergence in the total energy in order A 3 [because the inverse linear divergence 
seen in Eq. (|4.7b|) would cancel between the inside and the outside], and worse leads to a 
divergence in A 2 . It is probably not surprising that the situation is most satisfactory only in 
the conformally coupled case, because it is well known that any other choice leads to more 
singular quantum corrections |26|. 
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VI. CONCLUSIONS 



We have derived the Green's function in general for a massless scalar field for a spherically 
symmetric step-function potential. By taking the limit as the width of the step function 
goes to zero, we recover a ^-function shell potential, for which we consider both the weak 
and strong coupling limits. The latter corresponds to a Dirichlet shell, for which we recover 
the well-known results for the energy density and total energy. For weak coupling, we derive 
for the first time the behavior of the energy density as the surface is approached. We also 
examine the energy content of the shell itself, both for a thick shell, and in the thin shell 
limit. In the former case there are unphysical nonconformal divergences in the local energy 
density as the surfaces are approached which are exactly those found for a plane slab. These 
are most easily eliminated by using the conformal stress tensor. The inner and outer surfaces 
of a thick shell contribute a vanishing net surface energy, but there is a net effective surface 
energy in the thin shell limit, to be added to the integrated local energy density for the 
inside and outside regions, which is exactly the integrated local energy density of the shell. 

The situation is summarized as follows. The integrated local energy density inside, out- 
side, and within the shell is E in , E out , and E s , respectively. The total energy of a given region 
is the sum of the integrated local energy and the surface energy bounding that region: 

S m = £i„ + £-, (6.1a) 
£out = £out + <£+, (6.1b) 
£ s = E s + <B' + + <£'_, (6.1c) 

where (£± is the outside (inside) surface energy on the surface at r = a±, while <E'± is the 
inside (outside) surface energy on the same surfaces. Because for a nonsingular potential 
the surface energies cancel across each boundary, 

e + + C + = 0, <B + <£'_ = 0, (6.2) 

the total energy is 

£ = £ in + £ out + £ s = E in + E out + E s . (6.3) 

In the singular thin shell limit, the integrated local shell energy is the total surface energy 
of a thin Dirichlet shell: 

E s = + ^ 0. (6.4) 

This shell energy, for the conformally coupled theory, is finite in second order in the coupling 
(in at least two plausible regularization schemes), but diverges in third order. We show that 
the latter precisely corresponds to the known divergence of the total energy in this order. 
Thus we have established the suspected correspondence between surface divergences and 
divergences in the total energy, which has nothing to do with divergences in the local energy 
density as the surface is approached. This precise correspondence should enable us to absorb 
such global divergences in a renormalization of the surface energy, and should lead to further 
advances of our understanding of quantum vacuum effects. 
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